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Preface 



This book is based on the teaching documents used in the 
Physical Measurements Department at the University 
Institute of Technology in Saint-Nazaire, France, and is 
aimed at cooperative work/study engineering schools. It is not 
intended to provide a full grounding in signal processing for 
which readers should consult the reference works cited in the 
bibliography. 

I wish to thank my colleagues (M. Belkadi, T. Hairie, 
F. Poyer, H. Menana, J.-C. Olivier, B. Velay and Z. Shi) who 
have participated in this teaching activity and helped me to 
improve the quality of the content, and my students (in 
particular, A. Baudry, A. Bouvet, F. Gibon, R. Monnier and 
S. Vaillant) who, through their questions, have helped me to 
clarify and add certain exercises. 

The command files and programs mentioned in this book 
are available at http://www.univ-nantes.fr/auger-f (under 
"informations complementaires'Vadditional information). 



Frangois AUGER 
December 2013 



Introduction 



The aim of these practical exercises is to demonstrate 
some of the possibilities of practical implementation of signal 
processing techniques. To do this, we will use free software 
applications, including SoX, Audacity and Visual 
Analyser, focusing on their uses in real situations. 

SoX 1 is known as the "Swiss army knife of sound 
processing". Its main application lies in the handling of 
acoustic signals, but it may also be used for other types of 
signals, and nothing prevents its use for signals produced by 
any sensor (in addition to microphones). SoX is a free, which 
means that it may be used, studied and modified freely by 
users, and multiplatform, which means it may be executed 
using Windows, Linux or Mac OsX, software package. It is a 
command line application (see Figure 1.1), meaning that it 
does not include a human-machine interface using menus, 
buttons and fields, but is used by calling up the program 



1 Available for Windows, Linux and Mac OsX at http://sox.sourceforge.net. 
It was created in 1991 by Lance Norskog. The first version was known as 
AUX (for Aural Exchange). Since 1996, it has been maintained by Chris 
Bagwel. The program was in version 14.3.0 in January 2010, 14.3.1 in 
January 2011 and 14.4.0 in March 2012, showing the vitality and maturity 
of the development. 
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using parameters for the signals to process and the 
operations to apply. This method means that SoX may be 
used in application software written in any language. These 
command lines systematically take the general form: 



sox InputFile OutputFile Effect Ef f ectParameters 



Throughout these practical exercises, SoX may be used with 
Windows, either from a DOS command window 2 , or by 
creating and running command files 3 , using the .bat 
extension. These command files correspond to programs 
containing sequentially executed instructions. 



~ C:\WINDOWS\system32\cmd.ene 




D:SusersSf rancois\flctiuite EnseignenentSTpTsMp2010<DD>\SoxUorld>ren generation d 
e signaux elementaires 

D:\users\francoisSActiuite Enseignement\TpTsMp2010<DD>\SoxWorld>sox son3.mp3 -d 



son3 .np3 : 

File Size: 40.3k Bit Rate: 64.0k 

Encoding: MPEG audio 

Channels: 1 6 16-bit 

Samplerate: 48000Hz 

Replaygain: off 

Duration: 00:00:05.04 

In:99.Sy. 00:00:05.02 [00:00:00.02] Out:241k [ i ] Hd:0.0 Clip:?. 41k 

sox WARN sox: 'ossdsp' output clipped 7409 sanples; decrease volume? 
Done . 

D:\users\f rancoisSAct iuite EnseignementSTpTsMp2010<DD>\Soxlforld>pause 
Appuyez sur une touche pour continuer... 



Figure 1.1. Example of the use of SoX in a 

Windows command window 



2 In Windows, Marko Bozicovic's Console program, available at 
http://sourceforge.net/projects/console/, is useful. 

3 See http://en.wikipedia.org/wiki/Batchfile. For further information on the 
commands available using DOS/Windows, see http://en.wikipedia.org/wiki/ 
ListofMS-DOScommands. 
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The documentation for this program provides a succinct 
description of its possibilities. The aim of these practical 
exercises is to present certain signal generation, analysis and 
treatment functions offered by the program in order to 
correctly understand and use these possibilities. We will only 
consider processing operations that are not specific to the 
field of acoustics; readers who are interested in acoustic and 
musical effects may wish to carry out further, specific 
investigations. 

Ch is a multiplatform C/C++ interpreter designed for 
scientific calculations and for teaching algorithmics and 
computer programming. It conforms to ISO C 1990 and a 
large part of ISO C 1999. However, it only accepts part of 
C++. Ch was developed by Harry H. Cheng (University of 
California, Davis), who passed the responsibility of further 
development and commercialization to Softlntegration Inc. A 
"student" version is freely available to users who register 
with the website http://www.softintegration.com. 

Audacity 4 is another free, multiplatform software 
package. Its functions are sometimes comparable with, but 
generally complementary to, those offered by SoX: like SoX, 
it also offers generation, processing and analysis functions for 
acoustic signals. It operates in a graphical environment (see 
Figure 1.2), enabling easy visualization of signals. 



4 Available at http://audacity.sourceforge.net. 
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Figure 1.2. Main window of Audacity 



Finally, Visual Analyser is a free program 5 designed for 
Windows, which provides an oscilloscope, a spectrum 
analyzer, a low frequency (LF) generator and a frequency 
meter (see Figure 1.3). It is connected to the personal 
computer (PC) sound card by default, but may also be used 
with other acquisition cards. 



5 Available at http://www.sillanumsoft.com/. 



Introduction xiii 




Figure 1.3. Main window of Visual Analyser 
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Generating Elementary Signals 



1.1. General points 

A signal is a function of one or more variables, which has 
the specificity of being generated by a physical phenomenon 
and of carrying information rather than energy. Simple 
examples of signals include speech signals, acoustic signals 
and the electrical voltages provided by sensors. We may 
distinguish between different types of signals: 

- Continuous time signals, which have a (measurable) 
value at any time. 

- Discrete time signals, which only have (measurable) 
values at certain moments. 

- Sampled signals, which are time-discrete signals for 
which the time gap between two consecutive measurements 
is constant. 

- Deterministic signals, for which the value may be 
perfectly predicted using a mathematical model to calculate 
exact values at any time. 

- Random signals, which cannot be precisely predicted, 
either because it is impossible or because it is not necessary. 
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- Causal signals, which are null prior to the instant selected 
as time 0. These signals show the repercussions of a particular 
event occurring at time 0. 

- Stationary signals, which carry the same information 
throughout the observation period. All periodic signals are 
stationary, but not all stationary signals are periodic: for 
example, the deterministic signal defined by the expression 
x(t) = 10 sin(27r50t) + 5 sin(27r50v / 2t) is not periodic, as it is 
the sum of two sinusoids with frequencies in an irrational 
relationship. 

1.2. Generation of elementary wave forms 

To test the processing operations and analyses shown later 
in this book, it is useful to generate synthetic test signals 
(which are not the result of the observation of a real physical 
phenomenon). To do this, we use the synth command. The 
command file below (Generationl.bat) includes comment 
lines (which begin with the instruction rem) and six calls to 
the SoX program: 

rem generation de signaux elementaires avec SoX 
rem F. Auger, IUT Saint -Nazaire, dep. MP, dec. 2009 

sox -n sl.mp3 synth 3.5 sine 440 

sox -n s2.wav synth 90000s sine 660:1000 

sox -n s3.mp3 synth 1:20 triangle 440 

sox -n s4.mp3 synth 1:20 trapezium 440 

sox -V4 -n s5.mp3 synth 6 square 440 0 0 40 

sox -n s6.mp3 synth 5 noise 

for %%i in (200,300,400) do A 

sox -n s7_%%i.mp3 synth 15 sine %%i 

rem ecoute du resultat : le fichier de sortie est le canal 
rem de sortie par defaut (-d) , c'est a dire la carte son 

sox sl.rap3 -d 



pause 
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-The first command generates a file si .mp3 containing a 
recording of 3.5 s of a sinusoid with a frequency of 440 Hz. The 
-n (null) command is used to indicate the absence of an input 
signal. 

- The second command generates a file s2 . wav containing 
a recording of 90, 000 samples of a sinusoidal signal with a 
frequency that varies in a linear manner from 660 to 1,000 
Hz. 

- The third command generates a file s 3 . wav containing a 
recording of 1 min and 20 s of a triangular signal. 

- The fourth command generates a file s 4 . wav containing 
a recording of 1 min and 20 s of a trapezoidal signal. 

-The fifth command generates a file s5.wav containing 
a recording of 6 s of a square signal with a null offset (no 
continuous component), a null phase and a duty cycle of 50%. 
The command -V4 is used to choose the amount of information 
shown by SoX: its verbosity ranges from 0 (no information 
sent back to the user) to 4 (maximum information displayed). 

- The sixth command generates a file s6.wav containing a 
recording of 5 s of a white noise with a value ranging from 1 to 
1 with an average value of 0. 

- The seventh command uses a repetition structure to 
generate three files containing a recording of 15 s of a sinusoid. 
The variable %%i is used both for the file name and for the 
frequency of the sinusoid. The symbol A indicates that the 
instruction continues onto the following line. 

- The final command sends the contents of the file si .mp3 
to the default peripheral -d, i.e. the sound card, allowing us to 
listen to the signal. 

Exercise 1.1.- Write a command file to generate 10 files 
containing recordings of 2.5 s of sinusoids, for which the 
frequencies should be linearly spaced from 100 to 1,000 Hz. 
What is the analytical expression of these signals? 
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1.3. Elementary signal processing operations 

A number of elementary operations may be applied to 
existing signals. The command file below (called 

Generation2 .bat) shows some examples: 

rem transformations elementaires de signaux avec SoX 
rem F. Auger, IUT Saint -Nazaire, dep. MP, jan. 2010 

sox sl.mp3 sl_f aiblel .mp3 vol -6 dB 

sox sl.mp3 sl_f aible2 .mp3 vol -0.4 amplitude 

sox sl.mp3 s2.mp3 DeuxSonsl .mp3 

sox s2.mp3 -v 0.6 sl_f aiblel .mp3 DeuxSons2 . mp3 

sox sl.mp3 sl_avec_silencel .mp3 pad 1 

sox sl.mp3 sl_avec_silence2 .mp3 pad 1 0.5 

sox sl.mp3 sl_avec_silence3 .mp3 pad 1 5000s@3 0.5 

sox -m s3.mp3 s4.mp3 SommeSons . mp3 

sox sl.mp3 sl_avec_of f set .mp3 dcshift 0.05 

sox sl.rap3 sl_a_l_envers .mp3 reverse 

sox s3.mp3 s3_morceau . mp3 trim 1.5 2 

sox s3.mp3 s3_faded.mp3 fade t 10 1:00 20 

rem pause 

- The first two commands apply gains to the signal 
contained in file si .mp3; 6 dB in the first case (corresponding 
to a factor of 0.5) and 0.4 in the second case (corresponding to 
a reduction in amplitude of 60% and a phase shift of 180°). 

- The next two commands place two signals in sequence. In 
the second case, the amplitude of the first signal is reduced by 
a factor of 0.6. 

- The following three commands add null values (i.e. 
silence) to the signal contained in the file si .mp3. The first 
command adds 1 s of silence at the beginning, the second 
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adds 1 s at the beginning and 0.5 s at the end and the third 
additionally inserts 5,000 samples from the third second of the 
signal. 

-The following command mixes the two signals s3.mp3 
and s4 .mp3. 

- The next command adds a DC component equal to 0.05. 

- The following command reverses a signal (the end 
becomes the beginning, and vice versa). 

- The next command extracts a section of the signal 
contained in the file s 3 . mp3, beginning at 1.5 s with a duration 
of 2 s. This command is used to easily extract sections of 
signals for analysis or modification. 

- The final command is used to modify the amplitude of 
a signal. In this way, it is possible to vary the amplitude 
of a signal progressively, at the beginning and at the end. 
Parameter -t corresponds to a linear modulation. The second 
parameter indicates that the amplitude of the signal moves 
from 0 to its maximum value in the course of the first 10 s; the 
next parameter indicates that a decrease begins after 1 min, 
and the final parameter indicates that the linear decrease 
continues for more than 20 s. These last two parameters are 
optional. Other types of modulation are possible, as shown in 
the software user manual. 

Exercise 1.2.- Write a command file that generates a file 
we will call succession .mp3, containing the 10 signals 
generated in Exercise 1.1, placed in sequence. Next, write 
another command file to generate a second file, 
Sommesons .mp3, containing the sum of the 10 files generated 
in Exercise 1.1. What is the analytical expression of this 
second signal? What mathematical property does it satisfy? 

Exercise 1.3.- Write a command file to generate a 
recording of 1 min of signal equal to the sum of two sinusoids 
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with frequencies of 440 and 441 Hz and with the same 
amplitude. What does this signal look like? Justify its 
appearance using one of the following two relationships: 



Finally, write a command file to extract 3 s from the file 

AriaCantilenaVillaLobosPetibon .mp3 

beginning at 5 min 24 s, then listen to and visualize this 
extract using Audacity. How is this connected to the 
previous question? 

EXERCISE 1.4.- First generate a file Sinus .mp3 containing 
a sinusoid of frequency 440 Hz, then a file Noise. mp3 
containing random noise, both with a length of 7 s. Then, use 
a command such as 

sox -m Sinus. mp3 -v 0.01 Noise. mp3 SinusNoised.mp3 

to add the sinusoid to the noise, with the noise being 
attenuated by a factor of 0.01. Beyond what value of the 
attenuation factor is the noise no longer perceptible to the 
human ear? 

Finally, use a command such as 

sox -m Noise. mp3 -v 0.01 Sinus. mp3 SinusNoised.mp3 
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to add the noise to the sinusoid, with the sinusoid being 
attenuated by a factor of 0.01. Beyond what value of the 
attenuation factor is the sinusoid no longer perceptible to the 
human ear? 

Exercise 1.5.- The aim of this exercise is to show that it is 
possible to generate relatively complex files, which may then 
be used for tests or monitoring. 

Write a command file to generate an mp3 file containing a 
recording of 2 s of the sum of six sinusoids of frequency /o, 
2/o, 3/o, 4/o, 5/o and 6/0, where f 0 = 523 Hz. Next, generate 
two other files obtained using values of 659 and 830 Hz, 
respectively, for / 0 . Finally, generate a file containing the 
sequence of these three sounds, and then a last file that 
repeats the previous file three times 1 . 

For this exercise, a spreadsheet may be used with file 
GenerationCode.xls, which automatically generates a set 
of SoX instructions; these instructions may then be copied 
into a command file. 

Exercise 1.6.- Write a command file that: 

1) generates two sinusoidal signals with a duration of 10 
s and frequencies of 12,000 and 12,440 Hz. These two signals 
should be stored in files with extension . wav. 

2) calculates the sum of these signals and stores them in a 
third file with extension . wav. 



1 In 1968, Jean-Claude Risset used this principle to create a perpetually 
ascending sound, which gives the impression of constantly becoming higher 
and higher. 
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Listen to the obtained signal (preferably using headphones). 
What do you hear? Justify this result 2 . 

1.4. Using other file formats 

In the previous section, signals were stored in mp3 format 
files. However, SoX is able to handle a wide variety of file 
formats, both in terms of reading and writing. These formats, 
described in detail in the user manual, include . wav, used by 
the Windows operating system, and . dat, which corresponds 
to text files (coded in ASCII), in which the first line 
corresponds to the sampling frequency and the following lines 
contain two numbers for time and the value of the signal. The 
command file below (Generation3 .bat) 

rem generation d'un fichier texte au format .dat 
rem F. Auger, IUT Saint -Nazaire, dep. MP, jan. 2010 

sox sl.mp3 si. dat silence 0 trim 0.1 10s 

generates a file si . dat containing 10 successive values (10s) 
of file si .mp3 taken 0.1 s after the start of the signal (having 
removed the null values from the beginning of the file). The 
contents of file si . dat are as follows: 

; Sample Rate 4 8000 
; Channels 1 



0 



-0 . 68915808 
-0 . 65037082 
-0 . 60942747 
-0 .56646369 
-0.52161656 
-0 .47505816 
-0 . 42689275 
-0.37733861 



2.0833333e-005 
4.1666667e-005 



6.25e-005 



8.3333333e-005 
0.00010416667 



0.000125 
0.00014583333 



2 This technique is used by directional speakers to generate acoustic waves 
using ultrasonic transducers. See [MUI 11]. 
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0 . 00016666667 
0.0001875 



-0.32652545 
-0.27462651 



By creating files of this type, we are able to manipulate 
measurement signals of non-acoustic origin using SoX. We 
must simply normalize these signals so that they always 
fall between — 1 and 1. By choosing a sampling frequency 
suited to human audition, we are able to listen to 
non-acoustic phenomena. The program below (called 
GenerationFichierDat2010 . ch), written in C for the Ch 
environment, reads numbers from a text file, extracts the 
maximum and the minimum, and then generates a . dat file, 
applying cross-multiplication to the values of the text file to 
convert values between x m i n and x max into values between — 1 
and 1. The result may then be converted into an mp3 file, and 
may be listened to using the command 

sox -V4 Mesures.dat Mesures.mp3 



Programme permettant de transformer un fichier de mesures au format texte 
/ / en fichier de mesures au format . dat 
// F. Auger, janvier 2010 

#include <stdio.h> 

int main { ) 

{ 

char NomFichierTexte [ ] = "Mesures.txt" ; 

char NomFichierDat [ ] = "Mesures.dat" ; 

FILE ^PointeurFichierTexte, *PointeurFichierDat ; 

double x, xmax, xmin, xnorm, Somme, Difference, Fe=16000.0, Te~1.0/Fe, t; 

// on ouvre une premiere f ois le fichier pour trouver le min et le max 

PointeurFichierTexte = f open (NomFichierTexte, "r") ; 
if ( PointeurFichierTexte == NULL ) 
{ print f (" Impossible d'ouvrir le fichier %s en lecture \n", 



NomFichierTexte} ; } 



else 



// premiere lecture dans le fichier pour initialiser xmin et xmax 

if (fscanf (PointeurFichierTexte, "%lf ", &x) != EOF) 
{xmin-x ; xmax-x; } 
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while (fscanf (PointeurFichierTexte, "%lf ", &x) != EOF) 
{ 

if (x>xmax) {xmax=x; } 
if (x<xmin) {xmin=x; } 

} 

fclose (PointeurFichierTexte) ; 

Somme=xmax+xmiri ; Difference = xmax-xmin; t=0; 

// generation du fichier . dat par une relecture du fichier texte 

PointeurFichierTexte=f open (NomFichierTexte, "r" ) ; // ouverture en lecture 
PointeurFichierDat =f open (NomFichierDat , 
V);// ouverture en ecriture 

if { PointeurFichierDat == NULL ) 
{print f ( " Impossible d' ouvrir le fichier %s en ecriture\n", 
NomFichierDat) ; } 

else 

{ 

// on precise la frequence d' echantillonnage du signal de sortie, qui 
// n'est pas forcement la meme que celle du signal d' entree 

fprintf (PointeurFichierDat, " ; Sample Rate %8 . 2f \n" , Fe) ; 

// un seul signal = une seule voie = un seul canal 

fprintf (PointeurFichierDat , " ; Channels l\n" ) ; 

while (fscanf (PointeurFichierTexte, "%lf " , &x} ! - EOF) 

{ 

xnorm= (2*x-Somme) /Difference; 

fprintf (PointeurFichierDat, "%f %f \n" , t, xnorm) ; 
t=t+Te; 

} 

fclose (PointeurFichierTexte ) ; // c' est f ini, on f erme les fichier s 

fclose (PointeurFichierDat ) ; 

} 

} 

return 0 ; 

} 



Exercise 1.7.- If x is between x min and x max , which 
intervals must contain 



^min and 2 ^ *^min ^ rj 
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Provide a justification of the expression of xnorm used in 
the program. Next, modify the previous program to convert 
the file ecg.txt into . dat format. This file contains an 
electrocardiogram sampled at 720 Hz and acquired using a 12 
bit analog-digital converter 3 . This file may then be converted 
into .wav format using SoX, then visualized using 
Audacity. 




Figure 1.1. Result of the creation of an mp3 file based on electrical 
currents measured at the terminals of an asynchronous machine 



This format may also be used to process signals generated 
by SoX in other programs, for example in a spreadsheet 
program such as LibreOf f ice Calc 4 . 



3 Source: http://www.physionet.org/physiobank/database/aami-ecl3/. 

4 See http://fr.libreoffice.org/. 
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2.1. Measuring the temporal and frequency 
characteristics of a signal 

Three commands may be used to analyze a signal using 
certain specific characteristics: 

-The program soxi, supplied with SoX, provides 
information on the structure of a file. Thus, the command file 

(Analysel . bat) 

rem analyse d'un fichier avec soxi 

rem F. Auger, IUT Saint -Nazaire, dep. MP, jan. 2010 

soxi sonl . mp3 

soxi sonl.mp3 > sonl_info_soxi.txt 
pause 

provides the information shown below. With the first 
command, this information is shown on-screen. Using the 
second command, the information is stored in a file 

sonl_inf o_soxi . txt. 



Input File 
Channels 
Sample Rate 
Precision 
Duration 



' sonl .mp3' 
1 

48000 

16-bit 

00:00:03.55 



= 170496 samples 



~ 2 66.4 CDDA 
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sectors 
File Size : 28.4k 

Bit Rate : 64.0k 

Sample Encoding: MPEG audio (layer I, II or III) 

-The stat command in SoX determines certain 
characteristic parameters of a signal. Thus, the command 
file 1 (Analyse2 .bat) 

rem analyse d' un f ichier avec sox stat 

rem F. Auger, IUT Saint -Nazaire, dep. MP, jan. 2010 

sox sonl.mp3 -n stat 

sox sonl.mp3 -n stat 2> sonl_info_stat.txt 
pause 

produces the information shown below. With the first 
command, this information is shown on-screen. Using the 
second command, the information is stored in a file 

sonl_inf o_stat . txt. 



Samples 


read : 




169344 


Length 


(seconds) : 


3. 


,528000 


Scaled by: 


21474S 


33647.0 


Maximum 


amplitude : 


0. 


. 983058 


Minimum 


amplitude : 


-0 . 


, 977098 


Midline 


amplitude : 


0 . 


. 002980 


Mean 


norm: 


0 . 


,599965 


Mean 


amplitude : 


0. 


, 000013 


RMS 


amplitude : 


0 . 


, 669044 


Maximum 


delta : 


0. 


, 059710 


Minimum 


delta : 


0 . 


. 000000 


Mean 


delta : 


0 . 


. 034553 


RMS 


delta : 


0 . 


. 038529 


Rough 


frequency : 




439 



1 Warning: no space between "2" and ">". Specialists may wish to note that 
the normal output from SoX is stderr and not stdout. The syntax 2> 
allows us to redirect the stderr flow to a file, while the command > or 1> 
redirects the stdout flow. In the absence of redirection, the results appear 
on the screen. 
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Volume adjustment: 1.017 
Try: -t raw -U -1 

-Finally, the command stats in SoX provides other 
characteristic parameters of a signal. Thus, the command file 

(Analyse3 . bat) 

rem analyse d'un fichier avec sox stats 

rem F. Auger, IUT Saint -Nazaire, dep. MP, jan. 2010 

sox sonl.mp3 -n stats 

sox sonl.mp3 -n stats 2> sonl_info_stats.txt 
pause 

produces the information shown below. With the first 
command, this information is shown on-screen. Using the 
second command, the information is stored in a file 

sonl_inf o_stats . txt. 

DC offset 0.000013 
Min level -0.977098 
Max level 0.983058 
Pk lev dB -0.15 
RMS lev dB -3.4 9 

RMS Pk dB -3.44 
RMS Tr dB -3.89 
Crest factor 1.47 
Flat factor 0.00 
Pk count 2 
Bit-depth 29/29 
Num samples 169k 
Length s 3.52 8 

Scale max 1 . 000000 
Window s 0.050 



The SoX user manual provides information on the nature 
of these characteristic parameters, which are too numerous to 
discuss in detail here; some are also redundant. 
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EXERCISE 2.1.- Analyze and verify the results obtained 
using command file 2 below (Analyse20130 62 4 .bat), paying 
particular attention to the effective value of the signal. 

rem analyse de signaux elementaires avec SoX 

rem F. Auger, IUT Saint -Nazaire, dep. MP, june 2013 

sox -n ana_sl.wav synth 5 sine 440 

sox -n ana_s2.wav synth 5 triangle 440 

sox -n ana_s3.wav synth 4 square 440 0 0 40 

sox ana_sl.wav -n stat 2> AnaStat.txt 
sox ana_s2.wav -n stat 2>> AnaStat.txt 
sox ana_s3.wav -n stat 2>> AnaStat.txt 

rem pause 

EXERCISE 2.2.- What is the sampling frequency of the 
signals generated in Exercise 1.1? Write a command file to 
produce a highly detailed analysis of file SonRigolo3 .mp3. 
Discuss the obtained results. 

2.2. Time-frequency analysis of a signal 

The structure of certain signals may be clearly highlighted 
using a time-frequency representation. This representation 
corresponds to a distribution of the signal energy in a plane 
made up of the temporal and frequency dimensions. In this 
way, we are able to visualize the way in which the frequency 
representation of the signal evolves over time. The obtained 
image may be compared to a music score, with time shown on 
the x axis and frequency on the y axis. For example, the 
command file (Analyse4 .bat) 



2 The use of the command » allows us to place all results into the same 
file: the output flow is redirected to the end of the file AnaStat . txt. 



Signal Analysis 17 



rem analyse d' un f ichier avec sox spectrogram 

rem F. Auger, IUT Saint -Nazaire, dep. MP, jan. 2010 

sox -n s6a.wav synth 3 sine 660-2640 
sox -n s6b.wav synth 3 sine 1320-5280 
sox -n s6c.wav synth 3 sine 1980-7920 
sox -m s6a.wav s6b.wav s6c.wav s6.wav 

sox s6.wav -n spectrogram -o s6_sp.png 

sox s6.wav -n spectrogram -m -o s6_sp2.png 

sox s6.wav -n spectrogram -1 -o s6_sp3.png 

sox s6.wav -n spectrogram -1 -m -o s6_sp4.png 

sox s6.wav -n spectrogram -1 -m -S 0.5 -d 1.3 -o s6_sp5.png 

rem pause 



generates a file s6.mp3 containing a recording of 3 s of the 
sum of three sinusoids of which the frequencies increase over 
time. This signal is then analyzed using a spectrogram, which 
is a specific method of time-frequency analysis. The results of 
this analysis are stored in the output file (-o) s6_sp.png. 
The following command generates the same spectrogram in 
grayscale rather than in color, a format suitable for printing 
using a black-and-white printer. The next two commands use 
option -l, 3 which produces a white background. The final 
command analyzes the signal starting (-S) from instant 0.5 s 
over a duration (-d) of 1.3 s. The results of the penultimate 
command are presented in Figure 2.1. Among other things, 
this representation clearly shows the fundamental 
component and its harmonics, with frequencies that increase 
progressively over time. The SoX user manual gives a 
detailed description of the various parameters of the 
spectrogram command. 



3 Note that "1" is the lowercase letter L and not the number 1 (one). 
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Figure 2.1. Time-frequency representation of signal s6. wav, 
obtained using a spectrogram. dBFS stands for dB full scale and 
indicates that the different shades of gray correspond to the value in 
dB of the spectrogram, normalized by its maximum value 



Exercise 2.3.- Use spectrograms to differentiate the three 
signals generated by the following commands: 

sox -n chirpl.wav synth 3 sine 1000:20000 
sox -n chirp2.wav synth 3 sine 1000+20000 
sox -n chirp3.wav synth 3 sine 1000/20000 

Exercise 2.4.- Complete the command file from Exercise 
2.2 to construct a spectrogram of signal SonRigolo3 .mp3, 
both as a whole and for specific parts of the signal. Discuss 
your results. 

2.3. Non-parametric spectral analysis 

Spectral analysis, the representation of the distribution of 
energy of a signal in the frequency domain, is used in a wide 
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range of fields of application (acoustics, vibration mechanics, 
biomedical engineering, telecommunications, astronomy, 
earth sciences, etc.). Based on the Fourier transform, spectral 
analysis is an efficient tool for the analysis of stationary 
signals. In this section, we will consider its use for a number 
of signals to demonstrate the possibilities it offers, clearly 
highlighting the fact that the field of application of spectral 
analysis is far wider than the examples shown below. 

A possible starting point for the presentation of this 
technique is the Fourier transform of a signal x(t), which is 
the function of the real variable / (the frequency) with 
complex values defined by: 

/+oo 
x(t) e^ 2nft dt 
-oo 

As the signal is only observed and recorded over a finite 
length of time L from an initial instant to, we wish to 
approximate X(f) using the Fourier transform of the 
observed part of x(t), 

rto+L 

X(f) = / h(t - t 0 ) x(t) e- j2nft dt 
Jt 0 

= [ L h(u) x(t 0 + u) e~ ftnfu du 
Jo 

where h(t) is a weighting function chosen by the user and used 
to specify the importance of x(t) in the calculation of X(f). 
This integral is then calculated in an approximate manner by 
subdividing the interval [0; L] into N subintervals of the same 
width T s = L/N and approximating each integral using the 
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rectangle method: 



X(f) = 



,v_1 r(n+l)T s 

V / h{u)x{t 0 + u)e-^ fu du 

n=0 J » T » 



Te ^2 h(nT s ) x(tQ + nT s ) e 



- 3 2irfnT s 



n=0 



L N ~ 1 

— h(nT s )x(to + nT s )e 



-j2nfnT s 



n=0 



The most widespread practice is to calculate this 
frequency representation at frequencies / = mjf. This gives 
us the expression 



which is often calculated using an algorithm known as the fast 
Fourier transform (FFT), for which we graphically represent 
the square module \X(f)\ 2 (or its value in dB) as a function 
of the frequency The number of samples N thus determines 
the duration of the recording, L = NT S , and the discretization 
step of the frequency axis, A/ = F s /N. 

2.3.1. Analytical case studies 

Exercise 2.5.- We wish to calculate the Fourier transform 
of a recording of finite length of a sinusoidal signal with 
frequency fa, defined by x(t) = A cos(27r/ 0 1 + p). We presume 
that the recording of this signal is of length L, begins at 
to = —L/2 and ends at +L/2. A rectangular window will be 
selected as the "smoothing window" for the signal: 
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h(t) = 



1 if -L/2<t<+L/2 
0 else 



Show that the Fourier transform of this sinusoid segment 
may be written in the form: 



Specify the expression of H(f) and calculate H(0). If the 
amplitude A of the sinusoid is estimated by measuring H (q° , 
give a majorant of the absolute value of the relative error. This 
may be obtained using the following triangular inequality: 

\x\ — \y\ < \x + y\ < |x| + \y\ 

Discuss the practical applications of the obtained result. 

Exercise 2.6.- Let the signal x{t) be equal to zero at 
negative times and to e~ i//T at positive times: 



We wish to perform a spectral analysis on a recording of 
that signal, whose length is L. This entails calculating the 
Fourier transform of a signal y(t) equal to x(t) if 0 < t < L, 
and equal to zero otherwise. Calculate the Fourier transforms 



X{f) = \A e + w H(f -h) + \A e-*e H(f + / 0 ). 
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of x(t) and y(t). Write the relation which must be satisfied by 
the length L of the recording in order for the maximum 
relative error committed by confusing Y(f) and X(f) to be 
less than 1%. 

Exercise 2.7.- Calculate the Fourier transform of the signal 



where A e IR + and T 6 IR* + . From this, deduce the Fourier 
transform of the signal 



where A e H+, f 0 e iR+ and T e An example of this type 
of signal is shown in Figure 2.2. 

Exercise 2.8.- 

1) Recall the definition of the Fourier transform X(f) of a 
signal x(t). Show that ^(/) and ^p-(f) are, to within certain 
multiplication factors which we shall determine, the Fourier 
transforms of t x(t) and t 2 x(t). 

2) The Fourier transform of x(t) = e 2 is X(f) = 

\/2ire~ 27T f . By changing a variable, deduce from this the 

_ t 2 

expression of the Fourier transform of y(t) = e ^ , with T 
being a real, positive parameter characteristic of y(t). 




for t > 0 
for t < 0 ' 
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3 T=0.0125s and„*300Hz 

x10 0 




0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 



Time (s ) 

Figure 2.2. Signals x(t) (dotted lines) and y(t) (solid line) 
mentioned in Exercise 2. 7 

3) With reference to the previous two questions, calculate 
the Fourier transform of the Ricker wavelet (also known as 
the "Mexican hat wavelet" - see Figure 2.3), defined by 

h{t) = -jL- ( 1 - to I e" 2 ^ 
V ; V2^T \ T 2 J 

2.3.2. Practical application 

Exercise 2.9.- Using the command file below (which 
corresponds to AnalyseSpectralel .bat), we may generate 
a signal which is equal to the sum of two sinusoids with 
frequencies of 400 and 500 Hz. For how many samples must 
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the spectral analysis of this signal be calculated in order to 
distinguish between the lines of the two frequencies in the 
spectrum produced by Visual Analyser? Examine the 
effects of the choice of a weighting window and of the number 
of averages on the obtained representation. 

rem generation de signaux elementaires avec SoX 
rem F. Auger, IUT Saint -Nazaire, dep. MP, fev. 2010 

sox -n AnaSpecl .mp3 synth 3:00 sine 400 
sox -n AnaSpec2.mp3 synth 3:00 sine 500 

sox -m AnaSpecl. mp3 AnaSpec2.rap3 AnaSpec4.mp3 

rem pause 

°- 5 l ! ! ! ! ! ! ! ! ! 1 



0.4 - 



0.3 - 



0.2 - 



0.1 - 



0 



-0.1 - 




-0.2 1 1 1 1 1 1 1 1 1 1 

-5 -4 -3 -2 -1 0 1 2 3 4 5 

Time (s ) 



Figure 2.3. Graphic representation of the Ricker wavelet 
(with T = 1 s) studied in Exercise 2.8 
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EXERCISE 2.10.- Using SoX and Visual Analyser, 
analyze signal SignalStationnairel .mp3 and a second 
signal, SignalStationnaire2 .mp3, in order to identify 
the content as precisely as possible. Once you have specified 
the parameters of your spectral analysis, do not hesitate to 
use recording functions (settings/capture/capture 
spectrum), and particularly spectrum recording in a text file 
(File/Save Spectrum as text). 



3 
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3.1. Sampling 

A clear understanding of effective signal sampling 
conditions is essential for correct manipulation and 
processing of discrete-time signals. Sampling is an operation 
that converts a continuous-time signal x(t) into a 
time-discrete signal x[n], obtained by taking the value of x(t) 
at multiple instants over the sampling period: 
x[n] = x(nT/ Delta). It is therefore easy to pass from x(t) to 
x[n}. We will now consider the opposite situation, i.e. the 
reconstruction of x(t) based on x[n]. Using simple sinusoid 
signals, we can clearly demonstrate that sampling leads to a 
loss of information: 

- The signal x(t) = A cos(2tt(/ + F s )t + <p), where F s = 1/T S 
is known as the sampling rate, produces the following signal 
after sampling: 

x[n] = A cos(27r/T s n + 2irF s T s n + 92) 

= A cos(27r fT s n + ip) 

After sampling, this sinusoid of frequency f + F s will be 
considered as another sinusoid of frequency /. This applies 
to all sinusoids of frequency / + kF s , where A; is a positive 
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integer. This phenomenon, which is due to the periodicity of 
trigonometric functions, is known as frequency aliasing by 
translation. This phenomenon can cause the sampling process 
to modify the frequency of a sinusoid. To avoid this possibility, 
the frequency must be lower than F s . 

- The two sinusoid signals z±(t) = A cos (2-k (4f + A/) tj 

and z 2 {t) = A cos (2tv — A/) tj , corresponding to two 
sinusoids of the same amplitude and with symmetrical 
frequencies with respect to F s /2, produce the same discrete- 
time signal zi[n] = Z2[n] = (— 1)™ A cos(27rA/T s n) 
after sampling. This phenomenon is a result of the 
elementary properties of trigonometric functions (particularly 
the expressions of 003(6*1 + 9 2 ) and 003(6*1 — #2)), and is known 
as frequency aliasing by symmetry. To prevent modification of 
the frequency of a sinusoid as a result of this sampling process, 
the frequency must be lower than F s /2. 

On the basis of these two simple results, C. E. Shannon 
demonstrated that any given signal x(t) may be sampled 
correctly at sampling rate F s if the following two conditions 
are satisfied: 

-A frequency / max exists such that X(f), the Fourier 
transform of x(t), is null for any frequency / greater than / max . 
This means that the signal does not include energy located 
above /max, which is the highest frequency contained in the 
signal. 

- The selected sampling rate F s must be at least twice f ma , x : 
Fs — 2/ max . 

If these two conditions are verified, it is possible to calculate 
the value of the continuous-time signal x(t) at any instant t 
using the values of x[n] = x(nT s ), using an expression known 
as Shannon's interpolation formula. 
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In conclusion, for a signal to be correctly sampled at 
frequency F s , two conditions must be satisfied: 

-All of the information contained in the signal must be 
below F s /2. If this condition is not satisfied, we must either 
change the sampling rate or accept a loss of information. 

- No significant energy should be located above F s /2. If this 
condition is not satisfied, we must use a low-pass filter to 
verify Shannon's theorem. This is known as an anti-aliasing 
filter. 

3.1.1. Analytical case studies 

Exercise 3.1.- A microprocessor constantly repeats a cycle 
comprising a measurement of the voltage present at the input 
to an analog-to-digital converter (ADC) and a recording of the 
result in a digital-to-analog converter (DAC). Its operation 
rate is such that it can perform 2000 such cycles every 
second. 

A sinusoidal voltage of frequency 1500 Hz is then applied 
at the input to the ADC, using a low frequency (LF) signal 
generator. 

What is the frequency of the signal output by the DAC? 

Exercise 3.2.- Figure 3.1(a) shows the samples of a 
sinusoidal signal whose frequency is fo = 396 Hz. The 
time-period between two successive samples is 10 ms. Is this 
signal correctly sampled? Analyze the result obtained if 
need be. 

Figure 3.1(b) shows the samples of a sinusoidal signal 
whose frequency is fo = 112 Hz. This signal was recorded 
over 0.6 s, providing 601 points. Is this signal correctly 
sampled? Analyze the result obtained if need be. 
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Figure 3.1(c) shows the samples of a sinusoidal signal 
whose pulsation is lu 0 = 250 rad/s, sampled at a rate of 90 
points per second. Is this signal correctly sampled? Analyze 
the result obtained if need be. 




-10 1 1 1 1 1 1 1 1 1 ' 

5 10 15 20 25 30 35 40 45 50 

(b) 




10 20 30 40 50 60 70 80 90 100 

Time signals n 



Figure 3.1. Temporal representations of the discrete-time 
signals studied in Exercise 3.2 



EXERCISE 3.3.- By decomposing it into a Fourier series, it is 
possible to consider any periodical signal with period T 0 as 
the sum of a constant term (the continuous component) and 
sinusoidal components, of multiple frequencies with the base 
frequency fo = jr. If we sample at 300 Hz (without an 
anti-aliasing filter) a periodic signal whose base frequency is 
fo = 45 Hz, at what frequencies will its first 10 sinusoidal 
components be perceived after sampling? You can fill in the 
table shown in Figure 3.2, but justify your answers. What 
will the periodicity of this signal be after sampling? 

EXERCISE 3.4.- We wish to acquire a signal in the form 
x(t) = A sin(^) to measure its amplitude A. 



Signal Processing 31 



1) At what time can this parameter be measured directly? 

2) The signal is sampled with a sampling period T s , and 
we try to deduce A from the maximum value of all the x[n] 
values. Figure 3.3 shows the sampled signals obtained for 
different values of T s . In what circumstances will the error 
in the measurement of A be greatest, and why? 

3) Consider the case of the first column in Figure 3.3, i.e. 
where the two values closest to A are obtained at times \ — 
and \ + 1 f. We wish to find the value of T s for which the error 
committed in this case will be less than 0.5%. Deduce from this 
an inequality between T and T s , and another linking F e = jr 
and / = What do these inequalities become if we want the 
maximum error to be less than 0.1%? Draw a conclusion from 
the results obtained. 
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Figure 3.2. Table used in Exercise 3.3 



Exercise 3.5.- A sensor delivers a signal which contains 
pertinent data in a frequency interval [0,B]. This signal is 
first run through a low-pass anti-aliasing filter with cutoff 
frequency f c , and then sampled at a frequency F s . We choose 
f c = 1.5 B and F s = 3a B = 2a f c , where a > 1. The objective 
of this exercise is to determine the value of a for which the 
components of the signal situated in the interval [0, B] are 
least altered. 
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1) Let / be a frequency in the interval [0,1?]. What is 
the lowest frequency /' which will modify the frequential 
response of the signal with frequency / by superposition at 
that frequency by spectral aliasing? 

2) The anti-aliasing filter used is an iV-order Butterworth 
filter, the modulus of whose frequential response, 



yfl + (f/fc) 2N 

is a decreasing function on M+. What is the frequency / of the 
interval [0, B] which would be disturbed by the frequency /' 
which is least attenuated by the filter? 

3) At the frequency /' found in the previous question, we 
want the signal to be attenuated by a factor of 100. Deduce 
the value of a for TV ranging from 1 to 5. Draw a conclusion 
about the results obtained. 

4) The Tektronix 2642A analyzer uses a 6th-order 
Butterworth filter with F s = 2.56 f c . What is the attenuation 
at the frequency /'? 

Exercise 3.6.- For a digital oscilloscope, the sample rate 
is an important parameter. 1 In the technical specifications of 
the devices, it is always the maximum sample rate that is 
specified. However, in certain cases, the true sample rate may 
be considerably lower, for two reasons: 

- the number of acquisition channels. Often, digital 
oscilloscopes have only one ADC, which is therefore preceded 
by an analog multiplexer to take acquisitions on multiple 
channels. In twin-trace mode (two channels), the sample rate 
is therefore halved; 



1 See the article "Oscilloscopes: comment choisir la frequence 
d'echantillonnage?", Mesures, No 760, pp 38-40, December 2003. 
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- the size (or "depth") of the memory, i.e. the number of 
samples that the oscilloscope can record during an acquisition. 
The oscilloscope thus calculates the sample rate so that the 
whole of its memory is used for an acquisition. 



T =T/6 T =T/7 T =T/8 T =T/9 




0 0.5 0 0.5 0 0.5 0 0.5 



Figure 3.3. Sampled signals x[n] (*) obtained from the signal 
x(t) (solid line) in Exercise 3.4 for T s ranging from T/6 to 
T/17, with T = 0.5 s and A = 2.5 

Hence, for example, with a digital oscilloscope with a 
memory depth of 2 Mega Samples (2 x 1024 x 1024 = 2097152 
samples), 

- the recording of one signal with a timebase of 100 fi s/div 
will take 10 x 100 fis = 1 ms, because there are 10 horizontal 
graduations on the screen of the oscilloscope. The sample rate, 
i.e. the number of samples per second, will therefore be 

memory depth 
(timebase) x (number of graduations) 

2 x 1024 x 1024 

= = = 2.097 GHz 

10 _d 
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- the recording of two signals (twin-trace mode) with the 
same timebase will take the same length of time, but the need 
to store two signals in the memory instead of one, and the 
presence of the multiplexer before the ADC, mean that the 
sample rate is halved, so F s = 1.048 GHz. 

We want to use a digital oscilloscope, with a 24-Mega 
Samples memory depth and a maximum sample rate of 
20 GHz to visualize signals containing transitory phenomena, 
oscillating at 700 MHz. 

1) Is it possible to completely and correctly visualize this 
kind of signal on a single channel using a timebase of 
100 ^s/div? 

2) Is it possible to completely and correctly visualize these 
transitory phenomena if we take an acquisition on one channel 
with a timebase of 1 ms/div? 

3) Is it still possible if we take an acquisition on two 
channels with the same timebase? How will the transitory 
phenomena then appear? 

Exercise 3.7.- We wish to sample a real signal whose 
energy is located solely in a frequency band of bandwidth B, 
centered on the frequency f c , so between f c —B/2 and f c +B/2. 

1) Given the normal rules of sampling, how must we choose 
the sample rate F s so that this signal is correctly sampled? For 
f c = 10 kHz and B = 80 Hz, why is this solution undesirable? 

2) By application of the principle of frequency translation, 
what sample rate F s must we choose so that the component of 
the signal with frequency f c — B/2 appears in the same way as 
a constant signal (zero frequency)? 

3) If the sample rate used satisfies the relation obtained 
in the previous question, what additional relation must F s 
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and B satisfy so that there is no frequency overlap due to 
spectral aliasing? Why is this the optimum solution? Find a 
value of F s which is appropriate for the values of f c and B 
given previously. 

3.1.2. Practical applications 

Exercise 3.8.- The program below (which corresponds to 
file Echantillonnage . ch), written in C for the scientific 
computing environment Ch, generates a file in . dat format 
containing a 20 s recording of a sinusoid with a frequency 
that varies from ft = 400 Hz to // = 6, 000 Hz. 

// Programme permettant de generer la version echantillonnee 
// d' une sinusoide de frequence croissante 
// F. Auger, janvier 2010 

#include <stdio . h> 
#include <math . h> 

int main { } 

{ 

char NomFichierDat [ ] = "Echantillonnage . dat " ; 
FILE *PointeurFichierDat ; 

double MySignal, Fe, Te, t=0.0, Duree, fi, ff, Omegai, alpha ; 
double Amp=0.9, Phi, MyPi ; 

Duree=20 .0; // duree du signal en secondes 

Te=40e-6 ; Fe=1.0/Te; // periode et frequence d' echantillonnage 

fi=400.0; ff=6000.0; // frequences initiales et finales 

MyPi=M_PI; Omegai - 2*MyPi*fi; alpha-2 *MyPi* ( ff-fi ) /Duree ; 

/ / ouverture du f ichier en ecriture 

PointeurFichierDat = f open (NomFichierDat , " w" } ; 

if ( PointeurFichierDat == NULL ) 
{ print f ( " Impossible d' ouvrir le fichier %s en ecriture \n", 
NomFichierDat) ; 

} 

else 

{ 
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// on precise la frequence d' echantillonnage du signal de sortie 

fprintf (PointeurFichierDat, " ; Sample Rate %8.2f\n", Fe) ; 

// un seul signal = une seule voie = un seul canal 

fprintf (PointeurFichierDat, " ; Channels l\n H ); 

while (t<Duree) 
{ 

Phi=Omegai*t+0 . 5*alpha*t*t; // phase instantanee 
MySignal-Amp*cos (Phi ) ; // calcul du signal 

fprintf (PointeurFichierDat, "%f %f \n", t, MySignal); 
t=t+Te; 

} 

f close (PointeurFichierDat ) ; // c'est fini, on ferme 

} 

return 0 ; 

} 

The sinusoid signal takes the form x(t) = 0.9 cos(<^(t)), with 
tp(t) = uJit + at 2 /2. Calculate tl(t) = %(t) and How does 
the frequency vary as a function of time? Launch the program 
with a sampling period of T s = 40 fis and listen to the result. 
The . dat file may be transformed into a . wav file using the 
command 

rem conversion d'un fichier .dat en .mp3 

rem F. Auger, IUT Saint -Nazaire , dep. MP, jan. 2010 

sox Echantillonnage.dat Echantillonnage.wav 

Relaunch the program with a sampling period of 125 fis, 
and listen to the results again. Run a spectrogram analysis on 
the obtained signal. What has happened? How does this relate 
to Shannon's theorem? 

Exercise 3.9.- Use the program from Exercise 3.8 to 
generate a sinusoid signal with a duration of 20 s, for which 
the frequency increases linearly from 400 to 8, 000 Hz, 
sampled using a period of 170 fis. Convert the resulting file to 
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.wav format. Listen (using headphones) and produce a 
spectrogram of the results. Analyze and explain your results. 



3.2. Resampling 

SoX only handles sampled signals. All signals recorded 
using digital media are discrete-time signals. The program 
offers the possibility of resampling signals by changing the 
sampling rate; this constitutes one of the highlights of SoX. 
The operation is carried out using the rate command, as 
shown below: 

sox SonRigolo . mp3 SonRigolo2 . mp3 rate 32k 
soxi SonRigolo .mp3 > SonRigolo_info.txt 
soxi SonRigolo2 . mp3 >> SonRigolo_info.txt 



This produces a new signal SonRigolo2 .mp3 by 
resampling 2 the signal contained in file SonRigolo .mp3 at 
32 kHz, instead of the original sampling rate of 44.1 kHz. The 
two subsequent soxi commands produce the file below, 
which clearly shows the sample rate of the two files: 



Input File 
Channels 
Sample Rate 
Precision 
Duration 
File Size 
Bit Rate 
Sample Encoding 



' SonRigolo . mp3 ' 

1 

44100 
16-bit 
00:00:03.58 
28 . 6k 
64 . 0k 

MPEG audio (layer I, II or III 



157834 samples = 268.425 CDDA sectors 



Input File 
Channels 
Sample Rate 
Precision 
Duration 
File Size 



' SonRigolo2 .mp3' 
1 

32000 
16-bit 
00:00:03.60 
21 . 6k 



115200 samples ~ 270 CDDA sectors 



2 When using the MP3 standard, the authorized sample rates are 48 kHz, 
44.1 kHz, 32 kHz, 24 kHz, 22, 050 Hz, 16 kHz, 12 kHz, 11, 025 Hz and 8 kHz. 
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Bit Rate : 48.0k 

Sample Encoding: MPEG audio (layer I, II or III) 



Exercise 3.10.- What operations does SoX need to carry 
out in order to correctly resample a signal with an original 
sampling rate of 40 kHz using a new sampling rate of: (1) F s2 = 
10 kHz? (2) F s2 =16 kHz? 

Exercise 3.11.- 

— Analyze signal EmissionFrancelnter .mp3 to 
determine the number of samples it contains and the 
original sample rate. 

-Extract the first 2 min from this signal using the trim 
command, discussed in section 1.3, and save this extract in a 
. wav file. 

- Resample this new signal using smaller and smaller 
sampling rate values, and produce a table showing the number 
of samples and the size of the obtained file (in bytes) as a 
function of the new sampling rate. Analyze your results. At 
what sample rate does the signal become unintelligible? 

3.3. "Analog" filtering 

Stationary linear systems are useful for conditioning 
measurement signals, i.e. highlighting the information 
carried by the signals. Generally, these systems are as 
follows: 

-Dynamic: their evolution is governed by a differential 
equation (for continuous-time systems) or by a recurrence 
equation (for discrete-time systems). 

-Mono-variable: they have a single input value (known as 
excitation) and a single output value (the response). 
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- Causal: the value of the output at instant t depends solely 
on the input value at instants t' < t. This property translates 
the fact that effects cannot precede causes. 

-Linear: they verify the superposition principle, which 
states that the response of a system to a sum of elementary 
excitations is the sum of the responses of the system to each 
of these excitations taken individually. 

- Stationary: if the same excitation is given as input on 
two successive occasions, in strictly identical initial conditions, 
the two responses obtained will also be strictly identical. This 
property reflects the fact that the same causes produce the 
same effects, excluding the possibility of evolution, by wear 
and tear or aging, of the system. 

If a continuous-time system satisfies all of these properties, 
we may show that the output y(t) may be deduced from the 
input x(t) by an operation known as a convolution product: 



where h(t) is a characteristic function of the system known 
as the impulsional response. One important consequence of 
this result is that the response of a stationary linear system 
to a sinusoid is a sinusoid with the same frequency, amplified 
and phase shifted by terms that depend only on the system in 
question and the frequency of the sinusoid: 




h{r) x(t — r) dr 



if 



x(t) 



X eft* ft 



then 



y(t) 



H T {f)XeP*l\ 



with 



HAf) 



J^h(r)e-^dT 
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The proportionality coefficient between y(t) and x(t), 
Hjr(f), is the Fourier transform of h(t). 

More generally, if y(t) is the response to an excitation x(t) 
of a stationary linear system, then Yjr(f) = Hjr(f) Xjr(f), 
where Xjr(f) and Yjr(f) are the Fourier transforms of signals 
x(t) and y(t). This relationship demonstrates the possibility 
of conditioning measurement signals using linear systems, 
which behave, in the frequency domain, as a mask applied to 
the input signal in order to preserve information (taking 
Hjr(f) « 1) and to remove disturbances (using Hjr(f) m 0). 
These systems therefore allow us to separate information 
from noise, if their frequency locations are different. 

3.3.1. Analytical case studies 

Exercise 3.12.- An ideal low-pass filter, which leaves the 
signal totally unaffected in the passband, and completely 
eliminates everything outside of it, 



cannot be created, in reality With electronic components, it is 
only possible to create filters with a rational transfer 
function. Butterworth low-pass filters 3 are a group of 
technically feasible filters whose behavior is fairly close to 
that of ideal low-pass filters. The square modulus of their 
frequency response is equal to: 




l/l < fc 
1/1 > fc 



H(f)\ 2 



1 



l + (///c) 



3 Stephen Butterworth (1880-1958), an English engineer, first worked at 
the British National Physical Laboratory (NPL) and then at the Admiralty 
Research Laboratory until 1945, [BUT 30]. 
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This family therefore uses the following two values as 
parameters: 

- The cutoff frequency f c (f c £ M* + ). 

- The iV-order of the filter (TV £ ]N*). 

Hence, in order to solve the problem of designing a low-pass 
filter, we find the values of these two parameters that satisfy 
specifications compiled on the basis of the problem at hand, as 
follows: 

-For all frequencies between 0 and a given frequency /i, 
the frequency response must not have a square modulus less 
than a given value 1 — e\\ 

V/, |/| < /i => 1 - ei < \H(f)\ 2 < 1, where 0 < e x « 1 [3.1] 

This frequency f\ is the upper bound of the passband of 
the filter, and e\ is linked to the maximum deformation of the 
signal in that band. 

-For all frequencies above a given frequency / 2 , (/ 2 > /i), 
the filter's frequency response must have a square modulus 
less than a given value e 2 : 

V/, l/l > h 0 < ^(Z)! 2 < e 2 , where 0 < e 2 < 1 - ei [3.2] 

This frequency / 2 is the lower bound of the stopband of 
the filter, and e 2 is linked to the minimum attenuation of the 
signal in the stopband. 

The interval [/i,/ 2 ] is the transition band of the filter, 
within which range the designer imposes no constraints. 

The objective of this exercise is to show how to deduce the 
values of N and f c from the four values of /1, / 2 , ei and e 2 given 
by the technical specifications. 
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1) On the basis of equation [3.1], show that f c and N must 
satisfy the inequality: 



2) On the basis of equation [3.2], show that f c and N must 
satisfy the inequality: 



Which value of N should be chosen? How does the value of 



4) When the value of N is chosen, show that relations [3.3] 
and [3.4] necessitate that f c be in an interval [/ cm m, /cmax]- 
Show that the frequency V/cmin/cmax lies within this interval, 
and can therefore be chosen. What does this value become if 

ei = e 2 ? 

5) Deduce from this the parameters of a Butterworth filter, 
which satisfies the following technical specifications: 

/i = 30 Hz, ei = 0.01, h = 400 Hz, e 2 = 0.03 

What happens if f± is multiplied by 10? What happens if ei 
and €2 are divided by 10? 




[3.3] 




[3.4] 



3) Deduce from this the inequality: 




N evolve when j£ and/or e\ and e 2 vary? 



Exercise 3.13.- The objective of this exercise is to show how 
it is possible to create a bandpass filter with adjustable central 
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frequency, using two low-pass niters. This technique is known 
as bandpass filtering by quadrature modulation 4 . 

The principle is illustrated in Figure 3.5. Two signals x\(t) 
and x 2 (t), obtained by multiplication of the input signal x(t) 
by y/2 cos(oj 0 t) and y/2 s'm(u>ot) (where ujq = 27r/ 0 ), excite two 
identical filters with frequency response H(f). The outputs of 
these filters are multiplied by the same signals, so that yi{t) = 
zi(t) V2 cos(ojni) and y 2 {t) = z 2 (t) y/2 sinfwni), and then added 
to construct the output signal y(t) = yi{t) + y2{t). 

1) With the supposition that x(t) = X e- 727 ^*, calculate x\{t), 
zi(t) and yi(t). To do so, transform the term cos(uiot) into a sum 
of two complex exponentials. Verify that x±(t) is the sum of two 
complex exponentials of frequencies f + fo and f — fo, and that 
yi (t) is the sum of three complex exponentials. What are their 
frequencies? 



2) Calculate X2{t), z 2 (t) and y 2 {t) in the same way. Finally, 
deduce that y(t) = G(/)Ie^, where G(f) = H(f - f 0 ) + 
H(f + fo). 

3) Using the fact that a signal can be written in the form 
x (t) = f+™ X(f) df, deduce the expression of y{t). 

4) If H( f) is an ideal low-pass filter, 



where f c <C fo, summarily represent G(f). 

5) The two low-pass filters have real impulse responses iff 
V/ 6 M, H(f)* = H(—f). This is a necessary condition for 



4 For further details about this structure, see [HAZ 00]. 
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them to actually be feasible. Show that G(f) also has a real 
impulse response. 



a 




0 20 40 80 100 120 140 160 300 400 500 



Figure 3.4. Example of a problem of designing a low-pass filter: the 
filter (see Figure 3.4(a)) must not deform the components of the signal 
less than 30 Hz fej = 0.01, fx = 30 Hz), and must eliminate the 
components of the signal greater than 400 Hz (e 2 = 0.03, / 2 = 400 Hz). 
The solution chosen (N = 2, f c = 120 Hz, see Figure 3.4(c)), which 
reasonably closely approximates the ideal low-pass filter (dotted 
lines), does indeed satisfy the specifications in the passband (see 
Figure 3.4(b)) and in the stopband (see Figure 3.4(d)) 



3.3.2. Practical applications 

SoX allows easy filtering of signals using discrete-time 
equivalents of first- or second-order continuous -time filters. 
The command file below (Traitementl.bat) shows these 
possibilities: 

rem filtrage d'un signal avec les commandes sox 

rem lowpass, highpass, bandpass et bandreject 

rem F. Auger, IUT Saint -Nazaire, dep. MP, jan. 2010 
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sox SonRigolo3 .mp3 SonRigolo_lpl . mp3 lowpass -1 500 
sox SonRigolo3 .mp3 SonRigolo_hp2 . mp3 highpass -2 2000 0 . 8q 
sox SonRigolo3 .mp3 SonRigolo_bp2 . rap3 bandpass 200 5q 

sox SonRigolo3 .mp3 SonRigolo_br2 . mp3 bandreject 
200 0.6q 

sox SonRigolo_bp2 .mp3 SonRigolo_bp2b .mp3 gain -20 
rem pause 

-The first command applies a first-order (-1) lowpass 
filter to signal SonRigolo3 .mp3 with a limiting frequency 
of 500 Hz. The output signal of the filter is stored in 

SonRigolo_lpl . mp3. 

-The second command applies a second-order (-2) 
highpass filter to signal SonRigolo .mp3 with a limiting 
frequency of 2,000 Hz and a quality factor of Q = 0.8. 
The quality factor Q is equal to 1/(2 m), where m is the 
damping coefficient. The output signal of the filter is stored 
in SonRigolo_hp2 . mp3. 

-The third command applies a second-order bandpass 
filter to signal SonRigolo .mp3 with a central frequency of 
200 Hz and factor q = 5. The output signal of the filter is stored 

in SonRigolo_bp2 . mp3. 

- The fourth command applies a second-order bandreject 
filter to signal SonRigolo .mp3 with a central frequency of 
200 Hz and factor q = 0.6. The output signal of the filter is 
stored in SonRigolo_bp2 . mp3. 

- All of the filters presented above have a unitary maximum 
gain. To amplify or attenuate the signal, we may use the 
gain action in SoX. The fifth command applies a gain 
of —20 dB (i.e. an attenuation of factor 10) to the signal 
obtained using the bandpass filter, and stores the result in 
SonRigolo_bp2b .mp3. 
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Figure 3.5. Principle of a quadrature modulation filter, 
studied in Exercise 3.13. The rectangles containing x symbols 
are multipliers 



SoX also allows us to trace the modulus of the frequency 
response of these filters using the command -]-plot 
gnuplot. It generates a command file for the gnuplot curve 
tracing program 5 . SoX does not apply any processing 
operations to the input signal in this case, but the signal 
must still be identified in order for the program to obtain the 
sample rate. The command file shown below 
(Traitement2.bat), for example, produces the four 
frequency responses shown in Figure 3.6: 

rem trace de la reponse f requentielle des filtres obtenus avec 
rem les commandes sox lowpass, highpass, bandpass et bandreject 
rem F. Auger, IUT Saint -Nazaire, dep. MP, jan. 2010 



sox --plot gnuplot SonRigolo3 .mp3 

sox --plot gnuplot SonRigolo3 .mp3 

sox --plot gnuplot SonRigolo3 .mp3 

sox --plot gnuplot SonRigolo3 . mp3 



-n lowpass -1 500 > Courbel.plt 

-n highpass -2 2000 0 . 8q > Courbe2.plt 

-n bandpass 200 5q > Courbe3.plt 

-n bandreject 200 0 . 6q > Courbe4.plt 



pause 

The files generated by SoX, with extension .pit, are 
command files that enable gnuplot to produce graphics. By 



5 gnuplot is available at http://www.gnuplot.info. 
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default, gnuplot calculates 100 values of the represented 
function. If this proves to be insufficient for a satisfactory 
graphical representation of the function, an instruction of 
type set samples 300 may be used after the command: 

set ylabel 'Amplitude Response (dB) ' 

Unfortunately, the gnuplot command file generated by 
SoX does not allow us to zoom in on specific zones of the 
graphical representation. To do this, we need to open the 
gnuplot command file (for example Courbel.plt) using a 
text editor and remove the [-35:25] from command 

plot [f=10:Fs/2] [-35:25] 2 0 * logl 0 ( H ( f ) ) 

which imposes upper and lower limits on the y -axis . 



SoX :.w:ces i.n 1 '3 - «|.t:-;:v SCO :in-. «i=th:Hz: .0 :r;i-=-44 1 0C. SoX effect highpaSS gain-0 frequBnC/-20O0 O=0B (rate=44100) 




Figure 3.6. Frequency responses obtained using lowpass, 
highpass, bandpass and bandreject commands with the 
addition of the -plot gnuplot command 
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To obtain characteristic temporal responses for a filter, we 
may use and adapt the program below (TestFiltre.cn), 
which generates a signal containing an impulse, a step and a 
ramp, and then filter this signal using the appropriate SoX 
command. Figure 3.7 shows the signal produced by this 
program and the results obtained using command 

sox TestFiltre.wav TestFiltre2.wav lowpass -1 10 
// 

Programme permettant de generer un fichier pour tester des filtres 
// F. Auger, juin 2013 

finclude <stdio.h> 

int main ( ) 
{ 

char NomFichierDat [ ] = "TestFiltre.dat" ; 
FILE *PointeurFichierDat ; 
int i, NbPoints=2000; 

double Fe=16000.0, Te=1.0/Fe, t=0.0, Amp=0.9; 

// generation du fichier .dat par une relecture du fichier texte 

PointeurFichierDat =fopen (NomFichierDat, 
"w" ),"// ouverture en ecriture 



if ( PointeurFichierDat == NULL ) 
{print f (" Impossible d'ouvrir le fichier %s en ecriture\n", 
NomFichierDat) ; } 

else 

{ 

// on precise la frequence d' echantillonnage du signal 

fprintf (PointeurFichierDat , " ; Sample Rate %8.2f\n", Fe) ; 

// un seul signal = une seule voie = un seul canal 

fprintf (PointeurFichierDat , " ; Channels l\n"); 

for (i=0; i<NbPoints; i=i+l) 
{ fprintf (PointeurFichierDat, "%f %f \n", t, 0.0); t=t+Te; 

// 1' impulsion 

fprintf (PointeurFichierDat, "%f %f \n", t, Amp); t=t+Te; 

for (i=0; i<NbPoints; i=i+l) 
{ fprintf (PointeurFichierDat, "%f %f \n", t, 0.0); t=t+Te; 



// 1' echelon 

for (i=0; i<NbPoints; i=i+l) 
{ fprintf (PointeurFichierDat, "%f %f \n", t, Amp); t=t+Te;} 
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for (i=0; i<NbPoints; i=i+l) 
{fprintf (PointeurFichierDat, "%f %f \n", t, 0.0); t=t+Te; } 

// la rampe 

for (i=0; i<NbPoints; i=i+l) 
{ fprintf (PointeurFichierDat, "%f %f \n", 
t, Amp* (i+1) /NbPoints) ; t=t+Te;} 

for (i=0; i<NbPoints; i=i+l) 
{ fprintf (PointeurFichierDat, "%f %f \n", t, 0.0); t=t+Te; } 

fclose (PointeurFichierDat) ; 

} 

return 0 ; 

} 



Exercise 3.14.- Note that the transfer function of a second- 
order bandpass filter with a central pulse of lo 0 and damping z 
is written as: 



H{p) 



with ojq 
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The aim of this exercise is to study the variation of the 
actual central frequency and the bandwidth A/ of a 
second-order bandpass filter as a function of its quality factor 
Q and the desired central frequency fo- To do this, we trace 
the frequency response of the bandpass filters obtained using 
different values of the quality factor Q (for example Q = 2 n / 2 , 
for n from -2 to 3) and of fo (for example fo = 300 2 k for k from 
0 to 5), and then measure the actual central frequency and 
the bandwidth Af for each response in order to compare 
them with theoretical values. To do this, we may use two 
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nested loops as shown in the command file below 

(passe_bande . bat): 

rem F. Auger, IUT Saint -Nazaire, dep. MP, mars 2012 

rem %%i et %%j sont deux variables 

rem le symbole A en fin de ligne indique 

rem que la commande continue sur la ligne suivante 

for %%f in (300,600,1200,2400,4800,9600) do " 
for %%Q in (0.5, 0.707, 1, 1.414, 2, 2.8) do * 

sox --plot gnuplot SonRigolo3 .mp3 -n bandpass 
%%f %%Qq > bode_%%f_%%Q.plt 

rem pause 




Figure 3.7. Signal obtained using program TestFiltre . ch and 
response to the signal from a first-order "analog" low-pass filter with 
a limiting frequency of 10 Hz 



Exercise 3.15.- Filter the first signal obtained in Exercise 
1.2 using a band-reject filter with a central frequency of 
500 Hz and a rejection band of 100 Hz. Measure the 
amplitude of the successive sinusoids to deduce experimental 
values for the rejection bandwidth and the modulus of the 
frequency response of the filter, which should be added to the 
theoretical curve. 

Exercise 3.16.- Command 

sox SonRigolo .mp3 SonRigolo_lp2 . mp3 lowpass -2 500 3q 



may be used to apply the time-discrete equivalent of a second- 
order analog low-pass filter, with a static gain of 1, a natural 
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frequency of fo = 500 Hz and quality factor Q = 3 to signal 
SonRigolo .mp3. The output signal of this filter is stored in 
file SonRigolo_lp2 .mp3. Note that the transfer function of a 
second-order low-pass filter with a natural frequency of fo and 
a quality factor of Q is written as: 



where uo = 2ir fo and Q = i.e. z = jq. The modulus curve 
of the frequency response of this filter (expressed in decibels) 
may be obtained using the command 6 



sox --plot gnuplot SonRigolo. mp3 -n lowpass -2 500 3q > lowpass2.plt 



This generates a file that may be used with gnuplot to 
obtain a plot of the type shown in Figure 3.8. 

If Q is greater than V2/2, this filter presents resonance 
(i.e. a maximum higher than 1 has a non-null frequency) of 
frequency: 



6 The curve obtained with this command is constructed using 100 values 
from the frequency response. If this is not sufficient to obtain a satisfactory 
graphical representation, we may modify the .pit file by adding an 
instruction of type set samples 30 0 after the command set ylabel 
'Amplitude Response (dB)'. This instruction allows us to produce a 
curve using 300 values from the frequency response. Note that in order 
to zoom in on the graphic, we must simply remove " [-35:25] " from the 
instruction plot [f=10:Fs/2] [-35:25] 2 0 * logl 0 (H ( f ) ) . 



H(p) 
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SoX effect: lowpass gain=0 frequency=500 Q=3 (rate=44100) 



B -10 




Frequency (Hz) 



Figure 3.8. Frequency response of a second-order low-pass filter 
obtained using the lowpass command followed by the -plot 
gnuplot command 



The frequency response modulus for this frequency is equal 



to: 



G r = \H(finf r )\ 



2Q 2 
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Note that G r is independent of fa. The limiting frequency 
at -3 dB is equal to: 
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Calculate H(jlo 0 ). What are the limits of f r , G r and f c 
when Q tends toward infinity? What are the limits of f r , G r 
and f c when Q tends toward V2/2? Write a command file to 
use SoX in order to plot the frequency response of 
second-order low-pass filters obtained with two values of Q (2 
and 20) and two values of f 0 (1000 and 10,000 Hz). Using these 
four curves, measure the resonance frequency f r , the 
resonance gain G r , the gain at frequency /o and the limiting 
frequency f c . Compare these results to the theoretical values 
and discuss your findings. 

3.4. Digital filtering 

3.4.1. Analytical case studies 

EXERCISE 3.17.- To eliminate the high-frequency 
components of a discrete-time signal x[n], we can envisage 
using averaging filters. Let yi[n] and 2/2 [n] be the signals 
defined by 

yi[ n \ = - {x[n] + x[n — 1] + x[n — 2}} 

2/2 N = ^ i x [ n ] + x[n - 1] + x[n - 2] + x[n - 3]} 

Si is the system with input x[n] and output yi[n] and 52 the 
system with input x[n] and output 2/2 [n]. 

1) Calculate the static gains of these two systems and 
analyze the results obtained. 

2) Calculate the Nyquist gains of these two systems and 
analyze the results obtained. 

3) Express 2/1 [n — 1] as a function of x[n — 1], x[n — 2] and 
x[n — 3], and yi\n — 1] as a function of x[n — 1], x[n — 2], x[n — 3] 
and x[n — 4]. Deduce the expressions of yi[n] — yi[n — 1] and 
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2/2 W — U2[n — Deduce from this the frequency responses (or 
complex gains, or transfer functions) H\ (A) and H 2 (X) for these 
two systems. 

4) Calculate -ffi(g) and #2(3) and analyze the results 
obtained. 

5) Using an expansion limited to the 3rd third order of 
the sine function sin(x) « x — x 3 /6, calculate an approximate 
value for the cutoff frequency of these two filters, and draw a 
conclusion from this. 

Exercise 3.18.- On the basis of a discrete-time signal x[n] 
obtained by sampling with the period T s of a signal x(t), we 
construct a second signal y[n] using the relation 7 : 



y[n] = ay[n — 1] H — (x[n] - x[n - 1]) , 



where a is a real parameter between —1 and 1. 

1) Does this relation define a linear system? Does it define 
a steady-state system? Justify your answers. 

2) Calculate the static gain and the Nyquist gain. Analyze 
the results obtained. 

3) For a = 1/3, calculate the response of the system with a 
unit step. We will carefully explain why, in this case, y[— 1] = 
0. Analyze the results obtained, and make the link with the 
previous question. 

4) Calculate the frequency response of this filter, and 
recalculate the results for Exercise 3.2. Using SoX, plot the 



7 Such a system is called a DC blocking filter. 
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frequency response of this filter for different values of a, and 
compare the results obtained. 

5) Using, for instance, the expansions limited to the first 
order of the exponential and the function 1/(1 + x), 

e x = 1 + x + x e(x) and = 1 — x + x e(x), 

1 + x 

deduce that the slope of this transfer function at the origin 
is proportional to ^T^aJ- What can we conclude from this? 
Deduce an approximation of the cutoff frequency for this filter. 
Is the approximation obtained a minorant or a majorant of the 
true value? 

6) Figures 3.9 and 3.10 show the modulus and phase of 
this frequency response for multiple values of a. Characterize 
these responses and find the results from the previous 
question. Justify the use of this filter to eliminate the 
continuing undesirable components from certain signals. 

7) This filter is excited by a sinusoidal filter of frequency A. 
Plot the filter's response for A = 1/3 and for different values 
of A. Perform a significant number of tests, and analyze the 
results obtained. Find the link with the previous questions. 

Exercise 3.19.- Consider the discrete-time system that, 
based on a sampled signal x[n], constructs a new signal y[n] 
defined by the relation: 

y[n] = bx[n] + bx[n — 1] + ax[n — 2] — ax[n — 4] 
— bx[n — 5] — bx[n — 6], 

where a and b are two real parameters. 

1) Calculate the static gain and the Nyquist gain of this 
system. Interpret the results obtained. 
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Figure 3.9. Modulus of the frequential response of the filter 
studied in Exercise 3. 18 

2) Using the superposition theorem, deduce the system's 
response to a signal that fluctuates slightly around a constant 
xo, alternately taking the values x 0 + u x and xo — u x . For 
example, we could take x[n] = xq + u x (— l) n . 

3) Show that this system's frequency response can take the 
form 

H(X) = 2j {a sin(2vrA) + b sin(4vrA) + b sin(6vrA)} e' j6nX 

4) Using the expansion limited to the third order of sin(x), 
calculate a and b so that H(X) « j27rAe~ j67rA . We then obtain 
a filter that gives an approximation of the derivative of the 
signal x[n], with a delay of three samples, so y[n] ~ T s ^ ((n — 
3)T S ). 

Exercise 3.20.- We are now interested in a discrete-time 
system that, on the basis of a sampled signal x[n], constructs 
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Figure 3.10. Phase of the frequency response of the filter 
studied in Exercise 3. 18 

a new signal y[n] by the relation 8 : 

y [n] = (x[n] + 2x[n- 1] +3x[n- 2] +Ax[n- 3] 
+3x{n - 4] + 2x[n - 5] + x[n - 6]) 



1) Calculate the first terms of this filter's impulse response. 
Does this filter have a finite or infinite impulse response? 

2) Calculate the static gain of this filter, i.e. the ratio 
between the output signal and the input signal when these 
two signals are constant: x[n] = x[n — 1] = x[n — 2] = X and 
y[n] = y[n — 1] = y[n — 2] = Y. Interpret the results obtained. 



8 See[CIR05]. 
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3) Calculate the Nyquist gain of this filter, i.e. the ratio 
between the output signal and the input signal when these 
two signals are alternate series: x[n] = (— l) n X = —x[n — 1] = 
x[n — 2] and y[n] = (— l) n Y = — y[n — 1] = y[n — 2]. Interpret 
the results obtained. 

4) Calculate the frequency response H(A) of this filter, i.e. 
the ratio between the output signal and the input signal when 
these two signals are of the form x[n] = X z n and y[n] = Y z n , 
where z = e j2wX . 

5) From this, deduce that this frequency response can be 
written in the form 



Using the relations of symmetrization of the type 1 — x 2 = 

x^ 1 (x — x^ 1 ), deduce that 



6) Using a second-order Taylor approximation around zero, 



*w = ^( 1+ *- ,+ *- ,+ '- , ) a = (ifi^)) 



2 




sm(K x) 



K 2 -l 



,2 



K sin (a;) 



1 - 



x 



6 



calculate an approximate value of the cutoff frequency for this 
filter. Analyze the quality of that approximation using Figure 
3.11. 
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Figure 3.11. Representation of the modulus (linear scale) of 
the frequency response of the filter studied in Exercise 3.20, as 
a function of the normalized frequency X = fT e = f/F s 

Exercise 3.21.- We are interested in a discrete-time 
system that, based on a sampled signal x[n], constructs a new 
signal y[n] by way of the relation 9 

y[n] = ax[n] + bx[n — 1] + 2 cx[n — 2] + bx[n — 3] + ax[n — 4] 

1) Calculate the static gain of this filter, i.e. the ratio 
between the output signal and the input signal when these 
two signals are constant: x[n] = x[n — 1] = x[n — 2] = X and 
y[n] = y[n — 1] = y[n — 2] =Y. Interpret the results obtained. 

2) Calculate the Nyquist gain of this filter, i.e. the ratio 
between the output signal and the input signal when these 



9 See [HAM 73]. 



60 Signal Processing with Free Software 



two signals are alternate series: x[n] = (— l) n X = —x[n — 1] = 
x[n — 2] and y[n] = (— l) n Y = —y[n — 1] = y[n — 2]. Interpret 
the results obtained. 

3) Calculate the frequency response H(X) of this filter, i.e. 
the ratio between the output signal and the input signal when 
these two signals are of the form x[n] = X z n and y[n] = Y z n , 
where z = e j2wX . 

4) Using symmetrization relations of the type 1 + x~ 2 = 
x^ 1 (x + deduce that this frequency response can be 
written in the form: 

H(\) = A(X)e-^ x 
where A(X) = 2a cos(4vrA) + 2b cos(2vrA) + 2c 

5) To calculate a, b and c, we impose that A(0) = 1 (static 
gain equal to 1), ^4(0. 5) = 0 (Nyquist gain equal to 0) and 

(0) = 0. Deduce the values of a, b and c. What type of filter 
is thus created, and what is its cutoff frequency? 

Exercise 3.22.- We are interested in a discrete-time system 
which, based on a sampled signal x[n], constructs a new signal 
y[n] by the relation 

y[n] = K (x[n] — 2ax[n — 1] + 2ax[n — 3] — x[n — 4]) 

1) Calculate the first terms of this filter's response to 
an impulse at time zero. Does this filter have a finite or 
infinite impulse response? Does this impulse response verify 
a property of symmetry? 

2) Calculate the static gain of this filter, i.e. the ratio 
between the output signal and the input signal when these 
two signals are constant: x[n] = x[n — 1] = x[n — 2] = X and 
y[n] = y[n — 1] = y[n — 2] =Y. Interpret the results obtained. 
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3) Calculate the Nyquist gain of this filter, i.e. the ratio 
between the output signal and the input signal when these 
two signals are alternate series: x[n] = (—l) n X = —x[n — 1] = 
x[n — 2] and y[n] = (— l) n Y = —y[n — 1] = y[n — 2]. Interpret 
the results obtained. 

4) Calculate the frequency response H(X) of this filter, i.e. 
the ratio between the output signal and the input signal when 
these two signals are of the form x[n] = X z n and y[n] = Y z n , 
where z = e j2nX . 

5) Using symmetrization relations of the type 1 + x~ 2 = 
x^ 1 (x + x~ l ), deduce that this frequency response can be 
written in the form: 

H{\) = j A{\) e~^ x 
where A(X) = AK sin(27rA) (cos(27rA) - a) 

6) To calculate a and K, we impose that ^4(Ao) = 1 and 
^(Ao) = 0. Thus, we obtain a bandpass filter with central 
frequency An. Deduce the values of a and K. 

7) Figure 3.12 shows the modulus of the frequency response 
obtained for An = 0.2. What is the passband of the filter 
obtained? 

EXERCISE 3.23.- A system generates a signal y[n] on the 
basis of another signal x[n] by a relation of the form: 

y[n] = ho x[n] + hi x[n — 1] + /12 x[n — 2] 
+hi x[n — 3] + ho x[n — 4] 

1) Which relation must be verified by the coefficients ho, 
hi and h 2 so that y[n] is always null when x[n] is a constant 
signal, of the type x[n] = xq? 
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2) Which additional relation must be satisfied by the 
coefficients h 0 , hi and h 2 so that y[n] is always null when x[n] 
is a linear signal, of the type x[n] = vo T s nl 

3) Which additional relation must be satisfied by the 
coefficients ho, hi and h 2 so that y[n] is always equal to 70 T 2 
when x[n] is a parabolic signal, of the type x[n] = n 2 l 

4) Which additional relation must be satisfied by the 
coefficients ho, hi and h 2 so that y[n] is equal to joT^n when 
x[n] is equal to x[n] = n 3 ? 

5) Deduce from this the values of ho, hi and h 2 . What, then, 
is the value of the signal y[n] when x[n] = xq+vq T s n+-^^- n 2 + 

hTj 39 

6 n • 
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Figure 3.12. Representation of the modulus (linear scale) of 
the frequency response of the filter studied in Exercise 3.22, as 
a function of the normalized frequency A = fT s = f/F s 
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3.4.2. Practical applications 

SoX also allows us to carry out discrete-time filtering 
using the commands fir and biquad. The command file 
below (Traitement3.bat) illustrates the possibilities of 
these two commands: 

rem filtrage a temps discret a l'aide des commandes sox fir et biquad 
rem F. Auger, IUT Saint -Nazaire, dep. MP, jan. 2010 

sox SonRigolo .mp3 SonRigolo_f irl . mp3 fir 0 . 1 0 .2 0.4 0 .3 

sox SonRigolo . mp3 SonRigolo_f ir2 . mp3 fir CoeffsFIR.txt 

sox SonRigolo . mp3 SonRigolo_biq.mp3 biquad 0.6 0.2 0.4 1 -1.5 0.6 

sox — plot gnuplot SonRigolo . mp3 -n fir 0.2 0.2 0.2 0.2 0.2 > CourbeFir . pit 

rem pause 

-The first command uses the signal x[n] contained in the 
file SonRigolo . mp3 to create a new signal y[n], obtained 
using a third-order finite impulse response (FIR) filter of 
equation: 

y[n) = 0.1 x[n] + 0.2 x[n - 1] + 0.4 x[n - 2] + 0.3 x[n - 3] 

-The second command uses a second fir filter for which 
the coefficients are contained in file CoeffsFIR.txt. In this 
file, the coefficients are separated by spaces or placed onto 
separate lines (return/enter key). Lines that start with the 
character # are comment lines. An example of content from 
the CoeffsFIR.txt file is shown below: 

# filtre RIF du 4eme ordre 

# F. Auger, jan 2010 
0.35 

0.62 
0.51 
0.62 
0.35 

- The third command uses the signal x[n] contained in the 
file SonRigolo .mp3 to create a new signal y[n], obtained 
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using a second-order infinite impulse response (IIP) filter of 
equation: 

y[n] = 0.6 x[n] +0.2x[n- 1} +0Ax[n-2] 

+1.5 y[n- 1] -0.6y[n-2] 

- The fourth command traces the modulus of the frequency 
response of a fourth-order FIR filter of equation: 

y[n] = - (x[n] + x[n - 1] + x[n - 2] + x\n - 3] + x[n - 4]) 
o 

The result is shown in Figure 3.13. 
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Figure 3.13. Frequency response of the FIR filter obtained using the 
fourth command from file Traitement3 .bat. 



Exercise 3.24.- Filter a signal made up of a sequence of 
sinusoids with suitable frequencies using an averaging filter 
of length 10, with equation: 



y[n\ 



1 

To 



x [n — k] 
\k=0 / 
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Deduce experimental measurements of the modulus of the 
frequency response of this filter at the successive frequencies 
of the signal and compare them to theoretical values. Verify, 
experimentally, that the frequency response of this filter is 
equal to zero at frequency F s /10, and determine, 
experimentally, its cutoff frequency 

EXERCISE 3.25.- The command sine allows us to filter a 
signal using a specific type of digital fir filters. The command 
file below (Partiel2010c .bat) illustrates its possible uses: 

rem Filtrage avec la commande sine de SoX 

rem F . Auger, IUT Saint -Nazaire, dep. MP, mars 2010 

sox InputFile . mp3 OutputFilel . mp3 sine -3500 

sox InputFile . mp3 OutputFile2 .mp3 sine -3500 -n 43 

sox InputFile . mp3 OutputFile3 . mp3 sine 4000 

sox — plot gnuplot InputFile . mp3 -n sine -12000 -n 151 
> FreqResp.plt 

rem pause 

- The first command applies a low-pass filter with a cutoff 
frequency (at —6 dB) of 3, 500 Hz to the signal contained in file 
InputFile .mp3. The resulting output signal is stored in file 

OutputFilel . mp3. 

- The second command operates in the same way, specifying 
that InputFile .mp3 should be filtered using a filter of order 
42 (the output of the filter is then calculated from 43 values of 
the input signal). If this is not specified by the user, as in the 
first case, SoX selects the order of the filter itself. 

- The third command applies a high-pass filter with a cutoff 
frequency (at -6 dB) of 4,000 Hz. 

- The final command plots the modulus of the frequency 
response of a low-pass filter of order 150 with a cutoff 
frequency (at —6 dB) of 12 kHz. 



66 Signal Processing with Free Software 



Design a signal including a sequence of sinusoids with 
frequencies ranging from 200 to 1,400 Hz, by steps of 50 Hz. 
Each sinusoid should have a duration of 1.5 s. Next, filter this 
signal using a low-pass filter of order 215 with a cutoff 
frequency (at —6 dB) of 800 Hz. In what form should the 
relationship between signal y[n], obtained as output from this 
filter, and signal x[n], the input, be written? Measure the 
amplitudes of the sinusoids obtained as input and output 
using this filter and deduce experimental values of the 
modulus of the frequency response. Compare these values to 
the theoretical curve plotted using the program. 

We wish to filter a signal in order to keep only those 
components located between 310 and 1,240 Hz, with the 
exception of components between 490 and 510 Hz. Propose a 
solution using filters of this type and apply it to the signal 
generated at the beginning of this exercise to check that we 
obtain the desired frequency response. 
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